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It was shown that tunneling current flowing through a system with Coulomb correlations leads to 
charge redistribution between the different localized states. Simple model consisting of two electron 
levels have been analyzed by means of Heisenberg motion equations taking into account all order 
correlations of electron filling numbers in localized states exactly. We consider various relations 
between Coulomb interaction and localized electron energies. Sudden jumps of electron density at 
each level in a certain range of applied bias have been found. We found that for some parameter 
range inverse occupation in the two-level system appeared due to Coulomb correlations. It was 
shown also that Coulomb correlations lead to appearance of negative tunneling conductivity at 
certain relation between the values of tunneling rates from the two electronic levels. 
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I. INTRODUCTION 

Non-equilibrium Coulomb correlations can drastically 
influence on the local charge distribution in the vicinity 
of impurity complexes in nanometer tunneling junctions. 
Coulomb interaction results in significant changes of 
each localized state electron filling numbers and current- 
voltage characteristics of impurity complexes. Adjusting 
parameters of a tunneling contact one can obtain nega- 
tive tunneling conductivity caused by Coulomb correla- 
tions in a certain range of applied bias. There are sev- 
eral experimental situations in which Coulomb interac- 
tion values are of the order of electron levels spacing or 
even strongly exceed this value. It usually takes place if 
the distance between several impurity atoms or surface 
defects is comparable with the lattice scale, so coupling 
between their electronic states can strongly exceeds the 
interaction of these localized states with continuous spec- 
trum. Another possible realization is a quantum dot or 
two small interacting quantum dots on the sample surface 
weakly connected with the bulk states. Such systems can 
be described by the model including several electron lev- 
els with Coulomb interaction between localized electrons. 
Electronic structure of such complexes can be tuned both 
by external electric field which changes the values of sin- 
gle particle levels and by electron correlations of localized 
electronic states. One can expect that Coulomb correla- 
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tions in non equilibrium situation result in spatial re- 
distribution of localized charges and possibility of local 
charge density manipulation governed by Coulomb cor- 
relations. In some sense these effects are similar to the 
"co-tunneling" observed in [l|, 0. Moreover Coulomb 
interaction of localized electrons can be responsible for 
inverse occupation of localized electron states and neg- 
ative local tunneling conductivity in a certain range of 
applied bias. These effects can be clearly seen if single 
electron levels have different spatial symmetry. 

Great attention was paid to electron transport through 
a single impurity or a dot in the Coulomb blockade and 
the Kondo W\ regimes. These effects have been stud- 
ied experimentally and are up till now under theoretical 
investigation But if tunneling coupling is not 

negligible the impurity charge is not a discrete value and 
one has to deal with impurity electron filling numbers 
(which now are continuous variables) determined from 
kinetic equations. 

Non-equilibrium effects and tunneling current spectra 
in the system of two weakly coupled impurities (when 
coupling between impurities is smaller than tunneling 
rates between energy levels and tunneling contact leads) 
in the presence of Coulomb interaction were described by 
self-consistent approach based on Keldysh diagram tech- 
nique in [T]|.[T2j. In the present work we consider the 
opposite case when Coulomb coupling between localized 
electron states strongly exceeds tunneling transfer rates. 

We suggest theoretical approach based on the Heisen- 
berg equations for localized states electron filling num- 
bers taking into account all order correlators of local elec- 
tron density [l^. Tunneling current in a two-level sys- 
tem of spinless fermions with infinite value of Coulomb 
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interaction has been investigated in But obtained 

results do not take into account any nontrivial pair cor- 
relations in the system for finite Coulomb correlations. 
If one is interested in kinetic properties for the applied 
bias range larger than the characteristic energy of corre- 
lations between localized and band electrons in the leads 
then Kondo effect is unimportant. In this case for the 
finite number of localized electron levels one can obtain 
closed system of equations for electron filling numbers 
and all their correlators. It allows to analyze the role 
of Coulomb correlations in charge redistribution and in 
formation main features of I-V characteristics. 



II. THE SUGGESTED MODEL 

We shall analyze tunneling through the two-level sys- 
tem with Coulomb interaction of localized electrons 
FigUl The model system can be described by the Hamil- 
tonian H. 

ia k(7 

+ '^'tki{ctcyCia +h.C.) +'^tpi{c+^Cia- + h.C.) (1) 
kia pio- 



left (sample) and right (tip) leads of tunneling contact re- 
spectively. ifc(p)- tunneling transfer amplitudes between 
continuous spectrum states and localized states with en- 
ergies Ei. Operators c^(^^)lck{p) correspond to electrons 
creation/annihilation in the continuous spectrum states 
k{p). riia = c^Cicr-two-level system electron filling num- 
bers, where operator Cia- destroys electron with spin a on 
the energy level e^. 

Tunneling current through the two-level system is writ- 
ten in terms of electron creation/annihilation operators 
as: 

I = ha = ^ Ikia = ^ nka = 
icr ka 

= J2 ^fc*(< 4<tC»^ > - < c+ Cfc^ >) (2) 

ki(7 

Let us consider h = 1 elsewhere, so motion equation 
for the electron operators product c^^Cicr can be written 
as: 

dc~^ c ' 

+ Uij{njcr + rij^cy) ■ cl^Cia - tki ■ {riia- - fk) + 

+ V^Cfco-C/cV + ^'=J4cTC^<T = (3) 

k' =ik 

where 
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FIG. 1: a). Energy diagram of two-level system and 
b). Schematic spatial diagram of experimental realization. 
Coulomb energy Uij correspond to the interaction between 
electrons on different energy levels. 

Indices k and p label continuous spectrum states in the 



fk = 4aCk<T (4) 

In order to get equation for the tunneling current one 
has to multiply equation (jS]) by combinations of electron 
filling numbers operators n^i^j^j.^ in the following way: 

(I - ni_cr)(l - ?12-(t)(1 - n2a)Ck^Cla = 

= {{tkl ■ (nia -fk)+Yl *k'l4aCk'a + tk24„Ci„) ■ 
k'^k 

■(1 - ni_CT)(I - n2-CT)(l - n2a)} ■ {si - SkY^ 

(5) 

= {[tkl ■ {n\a - fk) + Y VlcLcfeV + tk2ctaCla) ' 
k'^k 

•ni_cr(l - n2-o-)(l - n2cr)} ■ {ei - £k + C^ii}"^ 

(6) 

a 

a' k' 

■n2a' (1 - '^i-<t)(1 - n^-^' )} ■ {ei - eu + 

(7) 
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ni-an^^' (1 - n2_^' )4aCia = 

a' k 



(8) 



{{tkl ■ {ni^ - /fe) + XI ^k'l4a(^k'a + i/c24aCl'^) ' 
■n2-an2ail - "-1-a)} ' {^1 - £k+ ^Uu}^^ 



(9) 



'T-l-cr"2-(Tn2(TCj,^Clcr = 
{{tkl ■ (ni^ - /fc) + 51 VlCforCj^V + tk24a(^la) ' 

•'T'l-ff"-2-<T"2CT} • {£1 -£k + Uii + 2Ui2y^ 



(10) 



The relation nf^ ~ Uia- was used in these equations. 

Neglecting changes of electron spectrum and local den- 
sity of states in the tunneling contact leads caused by the 
tunneling current we shall uncouple conduction and lo- 
calized electron filling numbers. This means also that 
we neglect any correlation effects between localized and 
band electrons - like the Kondo effect. 

It is easy to check: 



(1 - ni_cr)(l - "2-(t)(1 - n2a) + 

+ni_o.(l - n2-(T)(l - n2a) + 
+ "2^' (1 - '^i-t) (1 - ^^2-a' ) + 

+ Y ni-an2^' (1 - n^_y ) + 

(7 

+n2-an2cril - ^T-l-cr) + ni^^n2-crn2a- = 1 (H) 

Thus summing up the right- and left-hand parts of 
equations ([Sl fTO)) we get an equation for < c^^Cia >, 
which gives us after summation over k an equation for 
the tunneling current through the two-level system. To- 
tal current is a sum of two contributions: 



Ikla = rfei{(nicr) - ((1 - ?ll-ff)(l - n2-cr)(l - n2(T))/fc(£l) - 

- (ni_cr(l - n2-a){l - n2a)) ' fk{ei + Uii) - 

- ('^2cr(l - n2-a){l - ni_cr)) ' fki^l + U12) - 

- (n2-cr(l - ?^2(t)(1 - f^l-cr)) ' fki^l + U12) - 

- (ni_crn2cr(l - n2-cr)) ' /fc(ei + Un + U12) - 

- (ni_crn2-cr(l - n2cr)) ' /fc(ei + Un + U12) - 

- ('^2<T?T-2-cr(l - ni-^)) • /fc(ei + 21/12) - 

- {ni^„n2-an2a) ■ fk{£i + Un + 2U12)} + 

+ tkltk2l^0kC2aCla + 

+ Y^^'^^^k'Aa^'k'a) ■ 
k'^k 

1-^ (1 - - ^2-cr)(l - n2a) ^ _^ 

£1 - £k 

_^ ^ ni_cr(l - n2-cy){l - n2a) ^ 



+ ( 
+ ( 



£1 + Uii - Ek 
Ea' "2^'(1 -"1-'t)(1 -'^2-^'), 



£1 + U12 - £k 



£i + Un + Ui2 



£k 



^ ei + 2Ui2-£k ^ ^ 



ni_g.n2-CTn2cr 
£1 + Un + 2C/i2 - Ek 



)} (13) 



In what follows we shall neglect terms tkitk2i^okC2aCi(T 



and terms proportional ''''^^''^^^f^''''' " in the expression 
(|13p as they correspond to the next order perturbation 
theory in the parameter Relaxation rates r^i-pj^ = 
TT ■ • h'o are determined by electron tunneling tran- 
sitions from the two-level system to the leads k (sample) 
and p (tip) continuum states. z/o/c(p)-continuous spec- 
trum density of states. The main equation for the current 
(|13P includes mean electron filling numbers , pair and 
triple correlators for the localized states, which have to 
be determined now. Equations for the total electron fill- 
ing numbers nio- n2cr on the levels 1 and 2 can be found 
from the conditions: 



dt 

dn2cr 

dt 



= 1 



^7)1 — 



= 1 



k2a -r J-p2a 



In2ry — 



(14) 



where tunneling current /po- can be easily obtained 
from Ika by changing indexes k p. 

Pair filling numbers correlators can be found in the 
following way: 



Ika — Ikla + Ik2<y 



(12) 



Where expression for the tunneling current Ik2a can be 
obtained by changing indexes 1 -f-?- 2 in the equation for 
the tunneling current Ikia, which has the following form: 



dt 



dt 



(15) 



Full expressions which determine the system of equa- 
tions for pair filling numbers correlators through the 
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higher order correlators in the stationary case have the 
form: 

{ ) = (Ffei + Tpi + Tfcj + Tpj) ■ 

- {Tk^fk{e, + U,j) + Tp,fp{e, + C/,,)) ■ 

- [Tk.Usj + U,,) + T„U{e, + [/„■)) • 

- {Tkifk{e^ + i7,i + Urj) + rpjj,(ei + C/,i + Uij)) ■ 
■ {ni-an^„'{l~n,^_^,)) - 

- {Vki^k{e^ + 2U,j) + T.p,,fp{e, + 2C/,,)) • 

- (rfe,/fe(e. + C/,, + 2U,j) + rpjp(e, + Uu + 2C/,,)) • 

- {Tkjfkisj + + + rp,/p(ej + + C/,,)) ■ 

- (rfe,/fc(£, + 2C/,,) + Tp,fp{e, + 2C/,,)) • 

- {Tk,fk{e, + + 2U,j) + Tp.fpie, + U,, + 2[/,,)) • 

(16) 

High order correlators can be found in the similar way: 



(nia-ni^a-nj^a) (Note that the system of equations p4)) - 
([T^ allows to analyze magnetic regime with riia- ^ Ui-a 
as well). So the system of equations for the pair corre- 
lators Kii = {niaUi^a), K22 = {n2an2-a) and K12 = 
{nian2a) after substitution of equation ([T8)) to (fT6)) has 
the form: 

an ai2 ai3 \ / Kn \ 
021 022 a23 X K12 = F (19) 
031 032 033 / \ -K^22 / 

with coefficients an: 



ail ^ 1 

ai2 = 2-nJ{ei + Uii) -n'({ei + Uii + Ui2) - 

- 2 ~ • n^ie2 + U22 + U12) ■ $1 

1 1 

ai3 = -nf (ei + 2 ■ [/12) • $1 (20) 

a2i = -nf (£2 + 2 ■ C/12) • $2 

022 = 2 ■ n^(e2 + t/22) - ^^^(e2 + C/22 + C/12) - 

- 2 • ^Ji • (£i + Uii + C/12) • $2 

1 2 

a23 = 1 (21) 



dt ' ^ ^ dt + 



(17) 



So expressions for high order correlations have the 
form: 



{nj^nj_^n^_^,) = {Ykj ■ fkiej + Ujj + 2Uij) ■ 
■ {{ni^^n-j^} + {ni^aUj^^)) + 
+ Tfei • fk{ei + 2Uij) ■ {njafij^c,) + 
+ ^Pr fpi^J+Un+2U,,)■ 
■ {{iii-an.j„) + (ni_cr%--o-)) + 
+ • fp{ei + 2C/y) • (njo-rij-o-)} • 
• {Ffc, -{3 + /,.(£, + 2C/y)- 
- /fe(e» + C/,, + 2C/y)} + 
+ Tp, •{3 + /p(e, + 2C/,j) - 



/p(£^ + C/„; + 2C/,;,)}}- 



(18) 



We consider here paramagnetic situation n^^- 

^z— (7i ij^ia'^ja) — {j^ia'^j—o') and (jlia'^i — cr'^ja) 



V2 



(n^(e2 + C/12) 



ri + r2 

nl{e2 + 2 • c/12) • (1 + 2 . ^2)) 



£132 = 1 



(nf(ei + C/i2) 



Ti + r 



ni (£i + C/ii + C/12) •(1 + 4-^2)) 

r2 



(n^ (£2 + C^i2) 



ni'(e2 + C/2 + C/i2)-(l+4-yli)) 



ri + r2 
4 



a33 



Tl ^ T 



(n{ (£i + C/i2) 



ri + r2 

nf(ei+2-C/i2)-(l + 2-Ai)) (22) 



where Ti = Tki + Fp^. If we introduce tunneling filling 
numbers in the absence of Coulomb interaction nf{ei), 
nT{e^ + U,,)■. 



Tk^fk{s)+TpJp{e) 



(23) 



then coefficients $i and Ai can be found as: 
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nf (e, + Uu) - nfjej + Uu + -U^j) 
3 + (e, + 2 • U,/) - nf{e^ + Uu + 2- U^j) 

nf (g, + Uu + 2- Ujj) 

3 + nf{e, + 2 ■ t/y) - nf{e, + U,, + 2 ■ U,j) 

\ ■ nf (e, + - \ ■ nf (e, + U,, + _ 
3 + nf (e, + 2 • C/y) - nj{e, + i7,, + 2 ■ 
^ ■ nf (gj + 2 ■ L/jj) + \ ■ nf (e, + t/^j + 2 ■ U^j) 
3 + nf (e, + 2 • [/„■) - nf (e, + C/,, + 2 • C/y) 

(24) 



/ (ei + ?7ii) ■ ni,^ ^ 

F = nl{e2 + U22) ■ n2a 

\ Trtr;'^^^'^^ + ^i2)"-2o- + p^^n^(e2 + Ui2)ni^ ^ 

(25) 

Pair correlators Kij can be expressed through n^j^ 
from equations (|19II25I) . Substituting the solution for 
higher order correlators, obtained from equation (jl6|) and 
(fTS]) to equation (fT4|) one can find < riio- > and finally 
the tunneling current. 

The determinant of the system (jl9[) can turn to zero or 
even becomes negative for some choice of the parameters 
and consequently electron filling numbers of the two-level 
system can get negative values at some ranges of applied 
bias voltage. Such invalid system behaviour is the result 
of our approximation because we neglected the interac- 
tion between the two localized electronic states due to the 
electron transitions to the continuous spectrum states in 
the leads and back. To improve the results it is neces- 
sary to take into account corrections which can be found 
using next order perturbation theory in the parameter 
j^, retaining the terms tkitk2VokC2aCia in equation 
In this case the final equations for ni„ have additional 
nonlinear terms and can be schematically written as: 

nia ■ [All + Mi"-L) + "2(7 ■ {A12 + M2'T-D = n^{£i) 
n2cr ■ {A22 + i'2n\„) + nia- ■ (^21 + J^i^L) = '^'^(£2) (26) 

Coefficients Aij , /i^ and Ui have rather simple but cum- 
bersome form and depend only on the tunneling filling 
numbers and parameters of the tunneling contact. 

In this paper we shall not regarded this case. 

III. MAIN RESULTS AND DISCUSSION 

The behaviour of and I-V characteristics strongly 
depends on the parameters of the tunneling system. At 
first let us analyze the situation in which tunneling rates 
from both localized states to the leads are approximately 
equal = tfe(p)2- Figures [2E demonstrate behaviour 

of filling numbers and tunneling current obtained from 



kinetic equations for the different values of Coulomb en- 
ergies Uij and various electron levels location relative to 
the sample Fermi level in symmetric T^i ^ Tpi and asym- 
metric Tki ^pii^ki 3> Fpi) tunneling contact taking 
into account all order correlators. The bias voltage in our 
calculations is applied to the sample. So if both levels are 
above(below) the Fermi level all the specific features of 
charge distribution and tunneling current characteristics 
can be observed at the negative(positive) values of eV. 

In the case of both energy levels situated above (Figl^l 
FiglS]) or below (FiglSj Fig|6]) the sample Fermi level we 
observe charge redistribution between electron levels of 
reentrant character. When applied bias increases two 
possibilities for charge accumulation for large values of 
Coulomb energies Uij are realized in turn. Charge can 
be localized on both electron levels equally ni = 712 or 
mostly accumulated on the lower energy level (rti < ri2). 
In Fig[3] there are two ranges of applied bias where upper 
level become empty ni = (^2 < eV < ei and £2 + ^^12 < 
eV < £1 -I- U12) for large values of the Coulomb energies. 
Decreasing of the Coulomb energies leads to the situation 
when charge is mostly accumulated on the lower energy 
level (Figint), but ni ^ 0. In the particular range of 
applied bias £2 < eV < £1 -I- U12 one can find that the 
charge is completely localized on the lower energy level: 
ni = 0. 

Taking into account all order correlators gives us an op- 
portunity to investigate tunneling through the two-level 
system in the case of small Coulomb energies Uij ~ £i{j)- 
Figure [5] demonstrates how filling numbers and tunneling 
current dependencies change due to decrease of Coulomb 
energies for the symmetric tunneling contact Tki = Tpi 
(asymmetric contacts show the same tendencies). We 
demonstrate the case of both electron levels localized 
above the sample Fermi level. 

If Coulomb interaction is of the order of single elec- 
tron energies, three ranges of applied bias appear, where 
inverse occupation takes place: rii > 712 (FiglSJa) (£2 -I- 
2C/i2 <eV <ei + Uu, £1 + 2Ui2 <eV <£i + Uu + U12 
and £i-ht/ii-h2t/i2 < eV < e2 + U22+2Ui2). Such situa- 
tion exists due to the condition that system configuration 
with two electrons on the upper level and one electron 
on the lower level has lower energy than configuration 
with one electron on the upper level and two electrons 
on the lower level for the parameters shown in Fig|5]3. 
Further decreasing of the Coulomb energies (FigO:) re- 
duces the effect of inverse occupation and finally local 
charge mostly accumulates on the lower energy level as 
it should be. 

We obtain that the effects of reentrant charge redis- 
tribution is more pronounced for asymmetric contact if 
tunneling rates to the sample are lager than tunneling 
rates to the tip. 

It is necessary to mention that without Coulomb in- 
teraction one can find filling numbers for both electron 
levels to be simple step functions which correspond to the 
tunneling filling numbers n'^{ei) shifted from each other 
on the value £1 ~ £2. 
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The effect of inverse occupation due to the Coulomb 
correlations is more pronounced in a system with elec- 
tron levels positioned on the opposite sides of the sample 
Fermi level. (FigUlFiglTl). Without Coulomb interaction, 
when r^(p)]^ — rfc(p')2j the difference of the two levels oc- 
cupation numbers (rti — n2 ~ rfcirj,2 — ^pi^k2) turns to 
zero. Taking into account Coulomb correlations of local- 
ized electrons in the two-level system results in inverse 
occupation of the two levels at the wide range of applied 
bias voltage (FigglFiglTl). 

In Fig.(|4ji,b) there are three ranges of applied bias 
where inverse occupation takes place (ei + Un < eV < 
£2 + 2C/12, ei -t- 2Ui2 < eV < £2 + U22 + U12 and 
El -I- Uii + U12 < eV). It is clearly evident (Fig|3l3,,b) 
that when applied bias doesn't exceed the value ei -I- U12 
the whole charge is localized on the lower energy level 
[ni = 0). With the increasing of applied bias inverse oc- 
cupation takes place and localized charge redistributes. 
The effect of inverse occupation strongly depends on re- 
lation between tunneling rates. It is most pronounced in 
asymmetric contact with more strong tunneling coupling 
to the lead k (sample). But we have not found inverse 
occupation if the two-level system strongly coupled with 
tunneling contact lead p (tip) (FigHt). In this case with 
the increasing of applied bias upper electron level charge 
increases but local charge still mostly accumulated on the 
lower electron level. 

Decreasing of the Coulomb energies results in disap- 
pearing of the inverse occupation (Figl7}3,c) and local 
charge mostly accumulates on the lower energy level. 
This clearly demonstrates the role of Coulomb interac- 
tion in described charge distribution effects. 

Tunneling current as a function of applied bias voltage 
for different level's positions is depicted in (Figl^J-FiglTJl- 
f) (tunneling current amplitudes are normalized on 2rfe). 
For all the values of the system parameters tunneling 
current dependence on applied bias has a step structure. 
Height and length of the steps depend on the param- 
eters of the tunneling contact (tunneling transfer rates 
and values of Coulomb energies). If both levels are sit- 
uated below the Fermi level fFig |3l6t i-f) upper electron 
level doesn't appear as a step in current-voltage char- 
acteristics but charge redistribution takes place due to 
Coulomb correlations. 

For approximately equal tunneling rates for both lo- 
calized levels current-voltage characteristics are mostly 
monotonous functions. But some new peculiarities ap- 
pear if tunneling rates are essentially different. In Fig |8l9l 
we show some results for the case tk{p)i ^ tk{p)2 ■ In this 
case an interplay between "single electron" nonequilib- 
rium occupation effects and Coulomb correlation effects 
exists and at certain bias charge redistribution is accom- 
panied by negative differential conductivity. 

The case of both energy levels situated above the sam- 
ple Fermi level is shown in Fig. ([8]). If the tunneling 
transfer rate from the sample to the lower energy level 
is the largest in the system and the tunneling trans- 
fer amplitude from the lower energy level to the tip is 



the lowest one (Figl8^,c), we see, that local charge in 
the system is mostly accumulated on the lower energy 
level. Vice versa if the tunneling transfer rate from the 
sample to the upper energy level is the largest one and 
from the upper energy level to the tip is the lowest in 
the system (Figl5jD,d), one can find that local charge is 
mainly accumulated on the upper energy level and con- 
sequently inverse occupation takes place. But due to the 
Coulomb interaction three ranges of applied bias exist 
where local charge is mostly localized on the lower en- 
ergy level £2 < eV < ei, £2 + U12 < eV < £1 + U12 and 
£2 + U22 + U12 <eV <£i + Uii + U12. 

Inverse occupation also takes place when energy levels 
are positioned on the opposite sites of the sample Fermi 
level (Figin^) or when both energy levels are situated 
below the Fermi level (FiglHt)). In any case Coulomb 
interaction modifies single electron occupation behavior, 
changing with applied bias normal occupation to inverse 
one or vice versa. 

If we look at Fig.([S^), we find several ranges of ap- 
plied bias where the charge is distributed differently. 
These intervals depend on Coulomb interaction values: 
the whole charge is accumulated on the lower energy level 
{ni — 0) for eV < £i + Ui2; inverse occupation exists (lo- 
cal charge is mostly localized on the upper energy level) 
forei-hC/12 <eV < £2+U22+Ui2 andei-HC/ii-hC/12 < eV; 
charge is equally accumulated on both electron levels 
m = n2 if £2 + U22 + U12 < eV < £1 + Uii + U12. 

If both energy levels are situated below the Fermi level 
(Fig|9l3) there are similar ranges of applied bias in which 
charge is distributed differently (equally for £1 < eV < 
£2 + U12, inverse occupation if £1 + U12 < eV < £2 + 
U22 + U12 and £1 + Uii + U12 < eV and so on). 

The appearance of negative conductivity regions is the 
most essential feature of the tunneling characteristics, de- 
picted in Figs.(I51;,d, [nt,d ). We want to stress once more 
that formation of negative conductivity is an interplay 
between non-equilibrium effects connected with the tun- 
neling current and Coulomb correlations. 



IV. CONCLUSION 

We investigated tunneling through the two-level sys- 
tem with strong Coulomb interaction between localized 
electrons taking into account all order correlators of lo- 
cal electron density. It was shown that charge redistri- 
bution between electron states is strongly governed by 
the Coulomb correlations and is of reentrant type. Elec- 
tron filling numbers dependence on applied bias becomes 
quite different from that for non interacting electrons. 
Existence of the charge redistribution effects means that 
adjusting the applied bias one can control spatial redis- 
tribution of localized charges. So wide possibilities for 
local charge accumulation and charge switching exist for 
such systems 

Besides this, at certain values of Coulomb interaction 
of localized electrons one can obtain correlation induced 
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inverse occupation of the two-level system in different 
ranges of applied bias. Inverse occupation is mostly pro- 
nounced in asymmetric contacts with different tunneling 
rates to the sample and to the lead, and when one energy 
level lies below the Fermi level and another one - above. 

Changing the parameters of the tunneling contact 
(tunneling rates of each level to the leads) we can ob- 
serve interplay between two mechanisms responsible for 
non-equilibrium occupation of each level: tunneling cur- 
rent induced inverse occupation of two-level system at 
particular ratio between tunneling rates (which exists in 
the absence of Coulomb interaction) and inverse occupa- 
tion connected only with Coulomb interaction of localized 



electrons. 

We revealed that for some parameter range system 
demonstrates negative tunneling conductivity in certain 
ranges of applied bias voltage. Negative tunneling con- 
ductivity is revealed in asymmetric case Tki ^ Tpi (Fig|8] 
and Figini) and is more pronounced if both energy levels 
are situated above the Fermi level. When energy levels 
are situated on the opposite sites of Fermi level negative 
tunneling conductivity is much weaker and when both of 
them are positioned below the Fermi level it is negligible. 
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FIG. 2: Two-level system filling numbers a).-c). and tunneling current d).-f). as a function of applied bias voltage in the case 
when both energy levels are situated above the sample Fermi level. Parameters ei = 0.6, £2 ~ 0.3, U12 ~ 1.0, Un = 1.4, 
U22 = 1-5 are the same for aU the figures. a),d).rfci = Ffea = 0-01, Fpi = Fp2 = 0-01; b),e).rfci = Ffea = 0-03, Fpi = F^a = 0.01; 
c),f).Ffci = Ffc2 = 0.01, Fpi = Fp2 = 0.03; 
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FIG. 3: Two-level system filling numbers a).-c). and tunneling current d).-f). as a function of applied bias voltage in the case 
when both energy levels are situated below the sample Fermi level. Parameters ei — —0.1, £2 = —0.3, U12 = 1.0, Un — 1.5, 
U22 = 1.6 are the same for all the figures. a),d).Ffei = Ffc2 = 0.01, Fpi = Fp2 = 0.01; b),e).rfci = Ffc2 = 0.03, Fpi = rp2 = 0.01; 
c),f).rfci = rfc2 = 0.01, Fpi = Fp2 = 0.03; 
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FIG. 4: Two-level system filling numbers a).-c). and tunneling current d).-f). as a function of applied bias voltage in the case 
when one energy level is situated above and another one below the sample Fermi level. Parameters ei = 0.2, £2 = —0.3, U12 = 
1.0, Uii = 1.4, U22 = 1.7 are the same for all the figures. a),d).rfci = Ffca = 0.01, Fpi = Fp2 = 0.01; b),e).rfci = F^a = 0.03, 
Fpi = Fp2 = 0.01; c),f).Ffei = F^s = 0.01, Fpi = Fp2 = 0.03; 



n1 n1 ^— n1 




FIG. 5: Two-level system filling numbers a).-c). and tunneling current d).-f). as a function of applied bias voltage in the 
case when both energy levels are situated above the sample Fermi level. Parameters ei — 0.6, £2 = 0.3, Ffci = Ffe2 = 0.01, 
Fpi = Fp2 = 0.01 are the same for all the figures. a),d).[7i2 = 0.7, Un = 1.0, U22 = 1.1; b),e).?7i2 = 0.35, Un = 0.5, U22 = 0.9; 
c),f).[/i2 = 0.35, Uii = 0.45, 1722 = 0.55 
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FIG. 6: Two-level system filling numbers a).-c). and tunneling current d).-f). as a function of applied bias voltage in the case 
when both energy levels are situated below the sample Fermi level. Parameters ei — —0.1, £2 ~ —0.3, Fki ~ ~ 0.01, 
Fpi = = 0.01 are the same for all the figures. a),d).l7i2 = 1.0, Un = 1.5, U22 = 1.6; b),e).l7i2 = 0.4, Un = 0.5, U22 = 0.65; 
c),f).(7i2 = 0.15, Uu = 0.25, U22 = 0.4 
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FIG. 7: Two-level system filling numbers a).-c). and tunneling current d).-f). as a function of applied bias voltage in the 
case when one energy level is situated above and another one below the sample Fermi level. Parameters ei = 0.2, £2 = —0.3, 
Ffei = Tk2 = 0.01, Fpi = = 0.01 are the same for aU the figures. a),d).i7i2 = 1.0, Un = 1.4, U22 = 1.7; b),e)./7i2 = 0.1, 
Un = 0.25, U22 = 0.8; c),f).[/i2 = 0.1, Un = 0.15, U22 = 0.25 
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FIG. 8: Two-level system filling numbers a).-b). and tunneling current c).-d). as a function of applied bias voltage in the 
case when both energy levels are situated above the sample Fermi level for different values of tunneling rates. Parameters 
ei — 0.6, 62 — 0.3, Ui2 — 1.0, Uii — 1.5, 1722 = 1.6 are the same for all the figures. a),c).rfci = 0,06, Fpi — 0,05, Ffe2 = 0, 15, 
Fp2 = 0, 005; b),d).Ffci = 0, 15, Fpi = 0, 005, F^a = 0, 06, Fp2 = 0, 05 
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FIG. 9: Two-level system filling numbers a).-b). and tunneling current c).-d). as a function of applied bias voltage in the case 
when one energy level is situated above and another one below a),c). and both energy levels are situated below b),d). the 
sample Fermi level for different values of tunneling rates. Parameters Ffei — 0, 15, Tpi = 0, 005, rfc2 = 0, 06, Tp2 — 0, 05 are 
the same for all the figures. a),c). ei = 0.2, £2 = -0.3, U12 = 1.0, Un = 1.4, U22 = 1.7; b),d).ei = -0.1, £2 = -0.3, U12 = 1.0, 
Uii = 1.5, U22 = 1.6 



